Theory of inplane magnetoresistance in two-dimensional massless Dirac fermion 

system 
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We present the theory of the inplane magnetoresistance in two-dimensional massless Dirac fermion 
systems including the Zeeman splitting and the electron-electron interaction effect on the Landau 
level broadening within a random phase approximation. With the decrease in temperature, we 
find a characteristic temperature dependence of the inplane magnetoresistance showing a minimum 
followed by an enhancement with a plateau. The theory is in good agreement with the experiment 
of the layered organic conductor a-(BEDT-TTF)2l3 under pressure. In-plane magnetoresistsnce of 
graphene is also discussed based on this theory. 
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I. INTRODUCTION 

Since the discovery of unconventional integer quantum 
Hall effect in graphene, [l], which is a single atomic 
sheet of graphite, massless Dirac fermions realized in con- 
densed matter systems have attracted much attention. 
Under magnetic field, a remarkable difference between 
conventional two-dimensional electron systems and two- 
dimensional Dirac fermion systems appears in the Lan- 
dau level structure. In conventional electrons, the Lan- 
dau level energies are equally spaced. Meanwhile the 
Landau level energies in Dirac fermions with the Fermi 
velocity v are given by 



E n = sgn(n) 



hv 



(i) 



where n = 0, ±1, ±2, ... and £b = y/h/eB is the magnetic 
length. Q For the case of Dirac fermions, the Landau 
levels are unevenly spaced. What makes a crucial dif- 
ference compared to the case of conventional electrons is 
the existence of the zero energy Landau level that plays a 
central role for the unconventional integer quantum Hall 
effect. 

Massless Dirac fermion systems are not restricted to a 
purely two-dimensional system. The layered organic con- 
ductor cv-(BEDT-TTF)2l3 under pressure shows remark- 
able physical properties associated with a Dirac fermion 
spectrum. Q Theoretically it has been predicted that this 
system is a massless Dirac fermion system 0, 0] where 
the Fermi energy is at the Dirac point and the Dirac cone 
is tilted. (8rlToj This massless Dirac fermion spectrum is 
supported by first principles calculations. fTTi [l2j Experi- 
mentally the observation of the negative interlayer mag- 
netoresistance [lH supports the massless Dirac fermion 
spectrum. Application of the magnetic field decreases the 
interlayer resistivity. This negative interlayer magnetore- 
sistance is consistent with the existence of the zero en- 
ergy Landau level. [T3| The interlayer resistance decreases 



* morinari@yukawa.kyoto-u.ac.jp 



in proportion to the inverse of the applied magnetic field. 
This magnetic field dependence arises from the zero en- 
ergy Landau level degeneracy. 

In this organic Dirac fermion system, an intriguing 
inplane magnetoresistance was observed. [5[ Under mag- 
netic field, the inplane resistivity decreases gradually as 
the temperature T is decreased for T > 100K. After 
reaching a broad minimum around 100K, the resistivity 
increases and then shows a narrow plateau region around 
several Kelvin. After that the resistivity increases again 
as the temperature is decreased further. 

In this paper, we present the theory of the inplane 
magnetoresistance in massless Dirac fermion systems in- 
cluding the Landau level broadening effect due to the 
Coulomb interaction between Dirac fermions and the 
Zeeman energy splitting. We compute the inplane longi- 
tudinal conductivity by the Kubo formula using the Lan- 
dau level wave functions for massless Dirac fermions. The 
Coulomb interaction effect on the Landau level broaden- 
ing is computed by the random phase approximation. 
The result is consistent with the inplane magnetoresis- 
tance observed in a-(BEDT-TTF) 2 I 3 .@ The theory is 
also applied to graphene. 



II. MODEL 

For the description of two-dimensional Dirac fermions 
in the x — y plane, we introduce two component spinor 
field operator ip a (x,y) where a = ± denotes the spin. In 
graphene and a-(BEDT-TTF)2l3, there are two Dirac 
points in the Brillouin zone. We assume that Dirac 
fermions are degenerate with respect to these valley de- 
grees of freedom. We do not consider inter-valley in- 
teraction and focus on the single-valley properties. The 
Hamiltonian is given by rl — Ho + Vc , where 

Ho = dx f dyipl (x, y) hv (k x a x + k v a y ^ ip a (x, y) , 

(2) 



2 



with fc T 



-id x ^y and a XjV the Pauli matrices. 



The 



term Vc describes the Coulomb interaction between 
Dirac fermions, Vc = (1/2) VqPqP-q, where V q = 

q 

2ire 2 / (4ireoe\q\) with e the dielecric constant. Through- 
out this paper, we assume that the Fermi energy is at the 
Dirac point. We do not include the effect of the Dirac 
cone tilt in a-(BEDT-TTF) 2 I 3 (Ref. @) because it turns 
out that tilt is unimportant for understanding the main 
features of the inplane magnetoresistance of a-(BEDT- 
TTF)2i3 as we shall see below. 

In a magnetic field, the kinetic energy of Dirac fermions 
is quantized into Landau levels, Eq. ([1]). Taking the Lan- 
dau gauge A = (0, Bx) 1 the Landau level wave functions 
are represented by <J> njfe (x, y) = exp (iky) 4>n,h (x) / \fL^, 
where L y is the system size in the y-direction and 



(x) 



C'n 



-i sgn n 




h\ n \ ( — 



l \n\-l 



kit 



T + «i 

*-B 



(3) 



with Co = 1 and C n — l/\/2 for and sgnn = 1(— 1) 

for n > 0(n < 0) and sgnn = for n = 0. 

Here h n (£) are the eigenstates of the harmonic 
oscillator Hamiltonian — <9|/2 + S, 2 /2, h n (£) = 

J ff„(C)exp(-e 2 /2)/(2" /2 7r 1 /4Vd), with H n (0 the 
Hermitc polynomial. 

In terms of the Landau level wave functions, the 
field operator ip (T (x,y) is represented by ip a (x,y) = 
Y^, & n ,k ( x jU) c n,k,<r- Using this form, we find that the 

n,k 

Fourier transform of the density operator p(x, y) = 



ni ,n2,/c,cr 
x Fn 1 ,n 2 (Q) C ri 1 ,k,cr Cn 2,k+q y .(j, 

where the function F ni1l2 (q) is defined bvfl5l 

Fni,n 2 (Q) = CniC„ 2 [j| ni | j |„ 2 | (q) 



(4) 



+ sgn(nin 2 ) J| ni |- 



-i(q)]- (5) 



For ri\ > ri2, the function J ni .„ 2 (q) has the following 
form 



J (a)-J ni1 f -^-^ t 

/ n 2 f 2 
v rn-12 I y s 

Xi n 2 



(6) 



and J n2 ,ni(q) = [4i,n 2 (- q)]*- Here £™(x) are the asso- 
ciated Laguerre polynomials. 



III. THE COULOMB INTERACTION EFFECT 
ON THE LANDAU LEVEL BROADENING 

Now we compute the Coulomb interaction effect on the 
scattering rate of Dirac fermions that leads to the Landau 
level broadening. As we shall show below the tempera- 
ture dependence of the Landau level broadening gives 
rise to a broad minimum in the inplane resistivity that 
appears around T = T m - m . (For the case of a-(BEDT- 
TTF) 2 I 3 , it has been reported^ that T min ~ 100K.) Al- 
though it is easy to include the Zeeman splitting in the 
calculation of the Landau level broadening, we present 
the calculation for the spinless case because the interac- 
tion effect plays an important role at high temperatures 
where many Dirac fermions are excited from the zero en- 
ergy Landau level while the Zeeman spin splitting effect 
is negligible. 

The single particle Matsubara Green's function for 
the Landau level with the index n is G n (iu) v ) 
1 \iui v — E n — E„(ia; 1/ )], where w„ = (2is+l)irkBT is the 
fermion Matsubara frequency. Within the random phase 
approximation, the self-energy E„(iw I/ ) is described by 



k R T 



xfn.n' (q) F n ,, n (-q) G n > (iu) v + ifl v >) ,(7) 



where 

Dq (*fi„) 



2^1 



(q)-F« 2l ni (-q) 



m ,?i2 
f (E ni ) - f (E n2 ) 
iyiis E ni ~r" E n2 



(8) 



The summation over the boson Matsubara frequency f2„/ 
in Eq. ([7]) is carried out by using the spectral representa- 
tion of Vq(iVL v ) = Vq/ [1 — VqDq(iQ v )]. Performing the 
analytic continuation iui u — > Lu+iS with 5 an infinitesimal 
number and after some algebra, we obtain 



S n (u + iS) = ~ 47r 3^ 2 J de J dqqlm[Vq(e + iS)} 



X ^2 F n,n> (<l)Fn>,n (~q 



B J-oo JO 

n(e) + f(E n ,) 



OJ + iS — En 



(9) 



The imaginary part of the self-energy, — Im £„ (uj + id) , 
leads to the Landau level broadening. Instead, we use 
an approximate form, = — ImS„ (E n + iS). We do 
not attempt to compute this quantity in a self-consistent 
manner. The Coulomb interaction plays an important 
role if there are large numbers of excited Dirac fermions. 
However, the number of the excited Dirac fermions is 
suppressed at temperatures less than the Landau level 
energy gap. In such a regime, we may treat the Coulomb 
interaction perturbatively. 

Figure[T] (a) shows for different Landau levels where 
we set [HI v / 2/Bhv/£ B = lOK/T" 1 / 2 and e = 300 



3 



that were estimated [I?} from the analysis of the in- 
terlayer magnetoresistance in a-(BEDT-TTF) 2 l3. (Note 
that at ambient pressure a large dielectric constant that 
is the same order of magnitude as our value has been 
reported. [18[ ) 



In 

the 



the numerical calculation, we used 
recursion formula for the function 
\jn\j(n + k)\x k/2 exp (-x/2) L k n (x) instead of the 
recursion formula for the associated Laguerre polyno- 
mials because L^x) and the factorials can be huge for 
Landau levels with \n\ 3> 1. The summation with respect 
to the Landau levels is taken from n = —50 to 50. At 
temperatures below ~ 10 K, remain constant. This 
behavior is understood from the energy gaps created 
by the Landau level structure: the Coulomb interaction 
plays an important role when there are excited Dirac 
fcrmions to higher Landau levels. In order to excite 
Dirac fermions to higher Landau levels, the temperature 
should be larger than the energy gap created by the 
Landau levels. Thus, the characteristic temperatures are 
determined from the energy gaps between the Landau 
levels. As shown in Fig. QJa), Tq behaves remarkably 
differently while the other T% (n ^ 0) behaves sim- 
ilarly. At low temperature below 30 K the effect of 
the electron-electron interaction is rapidly suppressed 
because of the large energy gap between the zero-energy 
Landau level and the |n| = 1 Landau level. Reflecting 
this fact, r„ decreases as we increase the magnetic field 
because the Landau level energy gaps increase. 
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FIG. 1. (Color online) (a) Temperature dependence of 
at B — 10T for different Landau levels with 8 = 0.1 for a- 
(BEDT-TTF) 2 I 3 . (b)Temperature dependence of r£ at B = 
10T for graphene. 



Figure [TJb) shows for graphene where we take e = 
2.5 for the dielectric constant [l9[ and \j2j Bhv j Ib = 
400K/T" 1 / 2 for the Landau level structure parameter. Q 
Although the temperature dependence of r„ is different 
from Fig. QJa) because of the parameter differences, it is 
common that the n — Landau level component behaves 
differently as compared with the n=0 Landau level. Since 
the Landau level energy gaps between the n = Landau 
level and the n = 1 Landau level for graphene is about 
1000 K, the value of Tq is negligible in the temperature 
range shown in Fig. [IJb). 



This result is consistent with the experiment [20( sug- 
gesting that the zero energy Landau level is quite sharp in 
shape compared with the other Landau levels. Compared 
to a-(BEDT-TTF) 2 I 3 , the almost temperature indepen- 
dent region extended until ~ 120 K. This is because the 
Landau level energy spacing in graphene is larger than 
that in a-(BEDT-TTF) 2 I 3 . As shown in Fig. QTb) the 
interaction effect on T„ is negligible for T < 100 K due 
to the large separation between the Landau levels. 



IV. INPLANE MAGNETORESISTANCE 

Now we compute the inplane longitudinal conductivity 
a xx using the Kubo formula, 21 



x 7 n, a 

(E - E n , a f + VI 



dE 



dE 



E\n\+l,o) 



r 2 

1 N- 



[E — E 



-\n\-l.o) 



F 2 

-\n\-l 



(10) 



where / is the Fermi distribution function and the 
Zeeman energy splitting is included as E n>a - = E n + 
g/iBcrB /2. Here \ib is the Bohr magneton and we set 
g = 2. The scattering rate is assumed to be T n = r + r^\ 
where To is associated with impurity scattering. In the 
following calculation we take Tq = 2K that was estimated 
from analysis of the interlayer magnetoresistance data 
fl3| at low temperatures, [la] To reduce the numerical 
computation time we use Padc approximants for the tern- 
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FIG. 2. (Color online) The inplane resistivity for different 
magnetic fields with To = 2 K. The normalization parameter 
po is taken as po = p^flOOK) at B — I0T to compare with 
the experiment in Ref. Q]- 

perature dependence of r„ . For Landau levels with n^O 
we used the same Pade approximant for rf because 
with n ^ behave similarly as shown in Fig. [Ha). 

Figure [2] shows the inplane resistivity, p xx = 1/ a xx for 
different magnetic fields. Note that <j xy = because 
the Fermi energy is at the Dirac point. Here we as- 
sume particle-hole symmetry so that the Fermi energy 
is fixed to the Dirac point even at finite temperatures. 
The minima appear around T m j n ~ 100K. These minima 
appear because of the onset of the Landau level splitting 
effect: The Landau levels with \n\ < 10 are well sepa- 
rated each other. But those separations are unimpor- 
tant for T ~ 100 K because of the temperature broaden- 
ing effect due to the derivative of the Fermi distribution 
function in Eq. (fTD]). For T > 100 K, Landau levels 
with 1 7i | < 10 are almost continuously distributed be- 
cause \E n+1 -E n \< r n+ i + T„. For T < 100 K, we find 
that \Eiq±i — Eio\ > Tioii + Tio from the temperature 
dependence of r„. So the Landau level splitting effect 
appears for T < 100K. We computed a xx without includ- 
ing T^, and confirmed that the temperature dependence 
of p xx for T > 100K mainly arises from the temperature 
dependence of Y%. 

The appearance of a minimum at a characteristic tem- 
perature T m i n in the inplane magnetoresistance suggests 
that T m in is a crossover temperature from the interaction 
dominant regime to the almost non-interacting regime: 
for T > Tmin, the Landau level broadening smears out 
the Landau level energy spectrum. In this regime, the 
Landau level spacing is unimportant, and the electron- 
electron interaction, which requires the excitations from 
one Landau level to higher Landau levels, plays an im- 
portant role. By contrast for T < T m j n , the Landau 
level broadening is less than the Landau level spacing. 
Thus, the excitations from one Landau level to higher- 
Landau levels are suppressed. The characteristic tem- 
perature T m ; n depends on e, v, and B. Although there is 
no simple analytical formula for T m j n , one can determine 
T m ; n from the inplane magnetoresistance measurement. 



The same analysis can be applied to the surface states of 
three dimensional topological insulators. [22l [23| 

With decreasing the temperature from ~ 100K the 
resistivity increases because the number of Landau lev- 
els contributing to a xx decreases. Below 10K a narrow 
plateau region appears. If we compute p xx omitting the 
Zeeman energy splitting, we have a peak instead of the 
plateau and p xx approaches a universal curve that is in- 
dependent of the magnetic field. The peak position is 
scaled by \J~B. So the presence of the plateau is associ- 
ated with the Landau level splitting between n = and 
n = ±1. Namely, including the Zeeman energy split- 
ting transforms the peak to the plateau. For T < 2Tq, 
p xx turns to increase again, and then p xx approaches a 
temperature independent value. We note that for a con- 
ventional parabolic dispersion case p xx monotonically in- 
creases with decreasing the temperature because the Lan- 
dau levels are equally spaced. 

All features stated above are consistent with the ex- 
periment [B| except for T < 2Tq. In the experiment, p xx 
does not approach a temperature independent value for 
T < IK but increases further with decreasing tempera- 
ture changing the slope at a characteristic temperature 
r exp . This behavior suggests that there is an another 
Landau level split ting probably associated with valley 
splitting. In Ref. |2J], a Kosterlitz-Thouless transition 
scenario was proposed. We will investigate this point 
further in a future publication. 

Now we comment on the tilt of the Dirac cone. In 
a-(BEDT-TTF)2l3, theoretical calculations suggest that 
the Dirac cone is tilted. In the presence of the tilt of the 
Dirac cone, the Landau level wave functions are deformed 
[Tol | that leads to anisotropy of the resistivity. However, 
the features of the inplane magnetoresistance are unaf- 
fected by the tilt. The temperature dependence of the 
inplane magnetoresistance is determined by the Landau 
level structure. Since the tilt of the Dirac cone just leads 
to a modification of the overall factor of the Landau level 
energies and does not affect the Landau level structure 
qualitatively, [T3l the tilt is unimportant for the temper- 
ature dependence of the inplane magnetoresistance. 

Using the theory, we are able to understand some re- 
sults about a xx in graphene. Figure [3] shows a xx for dif- 
ferent r at B = 10 T. We computed a xx for B > 10 T 
as well (not shown) and found similar behaviors. The 
results with To > 10 K are in good agreement with the 
experiment [25j for B < 8 T. Experimentally Tq is esti- 
mated [l^| as To ~ 30 K. For clean samples with To, we 
should observe a peak associated with the Zeeman split- 
ting around T = 2Tq — psB. For To = 5 K, the peak 
appears around 2Fo — HbB ~ 3 K as shown in Fig. [3] 
In the experiment reported in Ref. [25| , a xx decreases at 
low temperatures for B > 10 T. To understand this be- 
havior, we need to assume that a valley splitting occurs 
as discussed in the literature 12611. 



5 



e 2 




100 



T(K) 



splitting, and the Coulomb interaction effect between 
Dirac fermions. The Coulomb interaction plays an im- 
portant role at high temperatures where Dirac fermions 
are excited from the zero energy Landau level. We found 
that the n = Landau level behaves differently compared 
to the other Landau levels. The features observed in a- 
(BEDT-TTF^Lj are consistent with our result except for 
T < IK where a valley splitting may play an important 
role. This theory has also been applied to graphene. We 
have found a consistent behavior with an existing exper- 
imental data and have predicted the presence of a peak 
structure of conductivity in clean samples. 



FIG. 3. (Color online) The temperature dependence of a xx 
for graphene for different Fo at B — 10 T. 

V. CONCLUSION 

In conclusion, we have investigated the inplane resis- 
tivity of Dirac fermions under magnetic held. We have 
included the Landau level structure, the Zeeman energy 
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